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ABSTRACT 

We observed the Seyfert 1 galaxy NGC 3516 twice during the flight of Astro-2 
using the Hopkins Ultraviolet Telescope aboard the space shuttle Endeavour in March 
1995. Simultaneous X-ray observations were performed with ASCA. Our far-ultraviolet 
spectra cover the spectral range 820-1840 A with a resolution of 2-4 A. No significant 
variations were found between the two observations. The total spectrum shows a red 
continuum, f u ~ z/~ 189 , with an observed flux of 2.2 x 10 -14 erg cm~ 2 s _1 A -1 at 1450 
A, slightly above the historical mean. Intrinsic absorption in Lyman (3 is visible as well 
as absorption from O vi AA1032, 1038, N v AA1239, 1243, Si iv AA1394, 1403, and C iv 
AA1548, 1551. The UV absorption lines are far weaker than is usual for NGC 3516, and 
also lie closer to the emission line redshift rather than showing the blueshift typical 
of these lines when they are strong. The neutral hydrogen absorption, however, is 
blueshifted by 400 km s _1 relative to the systemic velocity, and it is opaque at the 
Lyman limit. The sharpness of the cutoff indicates a low effective Doppler parameter, 
b < 20 km s . For b = 10 km s^ 1 the derived intrinsic column is 3.5 x 10 17 cm -2 . 
As in NGC 4151, a single warm absorber cannot produce the strong absorption visible 
over the wide range of observed ionization states. Matching both the UV and X-ray 
absorption simultaneously requires absorbers spanning a range of 10 3 in both ionization 
parameter and column density. 



Subject headings: galaxies: active — galaxies: individual (NGC 3516) — galaxies: nuclei 
— galaxies: Seyfert — ultraviolet: galaxies 



1. Introduction 



Only 3-10% of Seyfert 1 galaxies show intrinsic UV absorption in the resonance lines of highly 
ionized elements ( |Ulrich 198£ ). Of these, NGC 3516 has shown the strongest and most variable 
absorption lines ( |Ulrich fc Boisson 1983 ; Voit, Shull, & Begelman 1987; |Walter et al. 1990 ; Kolman 
|et al. 1993 ; Koratkar et al. 1996| ). NGC 3516 is unusual for a Seyfert 1 in other respects as 



well. While as many as half of all Seyfert l's show absorption by ionized material intrinsic to the 
source, characterized as a "warm absorber" (Nandra & Pounds 1994), the only Seyfert 1 besides 
NGC 3516 with equivalent neutral hydrogen columns exceeding 5 x 10 22 cm" 2 is NGC 4151 ( Kolman 



et al. 1993; [Nandra & Pounds 1994fc |Yagoob, Warwick, fc Pounds 1989fc |Yaqoob et al. 1993; ) . In 



NGC 4151, these variable columns range from 1 to 12 x 10 cm (Yaqoob, Warwick, & Pound 



1989; Yaqoob et al. 1993| ), and a similar range of variation has been observed in NGC 3516 (Kolman 
et al. 1993| ; |Nandra fc Pounds 1994] |Kriss et al. 1996| ). Extended X -ray emission has been seen 



in the nuclei of Seyfert 2 galaxies flWilson et al. 1992 ; Weaver et al. 1995), but among Seyfert 
l's the only examples are NGC 4151 QElvis, Briel, & Henry 1983] [Morse et al.TMEj) and possibly 



NGC 3516 ( Morse et al. 1995 ). NGC 3516 is also one of the rare Seyfert l's with an extended 
narrow-line region (NLR) having a biconical morphology ( |Ulrich Pequignot 198"0 ; Pogge 1989| ; 
Miyaji, Wilson, fc Pefez-Fournon 1992| ; |Golev et al. 1995| ). A gain, NGC 4151 is the only similar 
counterpart ( [Evans et al. 1993 ). This rare combination of strong UV and X-ray absorption and 
extended narrow-line emission suggests that they may be related phenomena. 

Extended narrow-line emission in Seyfert galaxies is commonly associated with photoionization 
by a collimated source of radiation. Biconical morphologies are most often found in Seyfert 2 
galaxies ( Pogge 1989| ; Evans et al. 1994; Schmitt &: Kinney 1996| ), and they suggest that our line 
of sight to the central source of radiation is obscured. In the context of unified models of Seyfert 
galaxies (see the review by |Antonucci 1993 ), Seyfert l's present us with a direct line of sight to the 
broad emission-line region (BELR) and continuum source, whereas our line of sight in Seyfert 2's 
is blocked by a torus opaque from the mid-infrared to at least soft X-rays. If the torus collimates 
the ionizing radiation, then biconical morphologies should not be observed in Seyfert 1 galaxies. 
Given this line of reasoning, Evans et al. (1993) proposed that the UV-absorbing material on our 
line of sight in NGC 4151, but not in the torus proper, might collimate the ionizing radiation. 

Although the strength of the UV and X-ray absorption in both NGC 4151 and NGC 3516 
suggests that the two absorbing mechanisms are related, it is not clear how. Kolman et al.'s (1993) 
simultaneous X-ray and UV observations of NGC 3516 were inconclusive due to a lack of variability. 
Common UV and X-ray absorption at much weaker levels in some active galactic nuclei has been 
successfully modeled with a single warm absorber ( |Mathur et al. 1994 ; Mathur, Wilkes, fc Elvis] 



1995 ), but the wide range of ionization states of the UV absorber in NGC 4151 is not compatible 
with the simplest warm absorber models ( [Kriss et al. 1995 ). 



The similarities of NGC 3516 and NGC 4151 in their UV and X-ray absorption and in their 
biconical NLR's prompted us to explore the far-UV spectrum of NGC 3516 shortward of 1200 A 
using the Hopkins Ultraviolet Telescope (HUT). Our goal was to search for further evidence that 
UV and X-ray absorbing gas in active galactic nuclei (AGN) is related to the collimation mechanism 
for the ionizing radiation. To improve our understanding of the relationship between the UV and 
X-ray absorbing gas, we also performed simultaneous X-ray observations using the Japanese X-ray 
satellite ASCA. In this paper we present the far-UV spectrum obtained with HUT. A companion 
paper (Kriss et al. 1996) discusses the ASCA observations. 
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2. Observations 

We observed NGC 3516 on two occasions during the Astro-2 mission, once for 1518 s beginning 
at 5:35:34 UT on 1995 March 11, and again for 2200 s beginning at 1:32:03 on 1995 March 13. 
Both observations were through a 20" aperture during orbital night when airglow is at a minimum. 
The basic design of HUT is described by Davidsen et al. (1992). Briefly, a 0.9-m mirror collects 
light for a prime-focus, Rowland-circle spectrograph. A photon-counting detector sensitive from 
820-1840 A samples the dispersed spectrum at a resolution of 2-4 A with ~0.52 A per pixel. 
Improvements to HUT, its performance during the Astro-2 mission, and our basic data reduction 
procedures are described by Kruk et al. (1995). The raw data were reduced by subtracting dark 
counts, correcting for scattered geocoronal Lya emission and subtracting airglow. We then flux- 
calibrated the spectrum using the time-dependent inverse sensitivity curves developed from on-orbit 
observations and model atmospheres of white dwarfs. Statistical errors for each pixel are calculated 
from the raw count spectra assuming a Poisson distribution and are propagated through the data 
reduction process. As there was no evidence for variability in the UV data, the two separate 
observations are weighted by their exposure times and combined to form the mean flux-calibrated 
HUT spectrum of NGC 3516 shown in Figure |]. 

To model the spectrum of NGC 3516 and measure properties of the continuum, emission lines 
and the absorption lines, we use the IRAFQtask specf it flKriss 1994a| ). We fit the continuum with 
a power law in f\. The brightest broad emission lines (O vi A1034, Lya, and C iv A1549) are well 
described by power law profiles while single Gaussian components are adequate for the weaker broad 
lines. The power law profile has a functional form F\ oc (A/A ) ±ct , where a = In 2/(1 + FWHM/2c) 



(c.f. NGC 4151, Kriss et al. 1992). Additional narrow Gaussian cores are required for Lya and He n 
A1640. Single Gaussian profiles are used for all absorption lines other than the Lyman series. We 
allow extinction to vary freely following a Cardelli, Clayton, & Mathis (1989) curve with Ry = 3.1. 

To model hydrogen absorption in the Galaxy and in NGC 3516, we compute grids of transmis- 
sion functions including transitions up to n = 50. Using Voigt profiles of varying column density 
and Doppler parameter, we then convolve the transmission with the instrument resolution. Galac- 
tic neutral hydrogen is fixed at zero redshift at a column density of 3.35 x 10 20 cm -2 ( [Stark et 



al. 1992 ) with a Doppler parameter of b = 10 km s _1 . A sharp, redshifted Lyman edge is readily 
apparent in the NGC 3516 spectrum. Neutral hydrogen intrinsic to NGC 3516 is permitted to vary 
freely in column density, redshift, and Doppler parameter. 

Our best fit yields \ 2 /» = 1652/1582 for 1671 data points between 916 A and 1800 A 
(we omit a region from 1207-1222 A surrounding geocoronal Lya). The fitted continuum has 
fx = 3.61 x 10~ 14 (A/IOOOA) -0 - 11 erg cm' 2 s" 1 A" 1 with E(B-V) = 0.06 ±0.01. In frequency 
space this corresponds to a spectral index a = 1.89 for f v ~ v~ a . This is rather steep, but it is not 



1 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Observato- 
ries, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 
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Fig. 1. — The flux-calibrated spectrum of NGC 3516 obtained with the Hopkins Ultraviolet 
Telescope during the Astro-2 mission is shown. Significant emission and absorption lines are marked. 
The indicated Lyman limit is at a redshift of 0.0075, and it is intrinsic to NGC 3516 as are the 
absorption lines of O vi, N v, Si iv, and C iv. The remaining absorption lines are likely galactic 
in origin. The earth symbol indicates the residual feature produced by subtraction of geocoronal 
Lya emission. 



that unusual for a low-luminosity AGN. For comparison, the Astro- 1 spectrum of NGC 4151 showed 
a spectral index of 1.50 ( Kriss et al. 1992] ), the FOS spectrum of NGC 1566 has a = 1.53 flKriss et 



al. 199 1|) , and the UV continuum of M 81 has a = 1.5 - 2.0 ( [Ho, Filippenko, fc Sargent 1996|) . Our 
fitted extinction is higher than the predicted Galactic reddening of E(B — V) = 0.02-0.03 in the 
maps of Burstein & Heiles (1982), but compatible with that expected using Njji = 3.35 x 10 20 cm" 2 
and a gas-to-dust ratio of Nhi/E{B — V) = 5.2 x 10 21 cm" 2 ( [Shull fc Van Steenberg 1985|) . Previ- 
ous work based on IUE observations ( [Kolman et al. 1993 ; Koratkar et al. 1996|) have determined 
E(B — V) = 0.15 based on the strength of the 2200 A absorption feature, but such a large extinction 
correction provides a poor match to our data — for E(B - V) = 0.15, x 2 = 1782. To show the 
sensitivity of the data to the extinction correction, Figure ^ compares the observed spectrum with 
power-law continua after correction for E{B — V) = 0.06 and 0.15. The figure clearly shows that 
no extinction correction leads to a deficit in continuum flux below 1000 A, whereas an extinction 
correction of E(B — V) = 0.15 produces an excess in flux at short wavelengths. We suggest that 
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Fig. 2. — The figure shows the HUT spectrum of NGC 3516 on a log-log scale with extinction 
corrections applied as indicated. The spectra have been scaled up for display by the factors under 
each spectrum. The solid straight lines through each spectrum are the best fit power law continuum 
levels for the given extinction correction. No correction clearly produces a deficit in flux at short 
wavelengths, and E(B — V) = 0.15 leads to an excess in continuum flux shortward of 1000 A. 
E(B — V) = 0.06 provides the best fit. Note that the broad wings on the emission lines leave clear 
continuum visible only in narrow windows shortward of 980 A and near 1130 A, 1450 A and 1800 

A. 



the apparent strength of the 2200 A dip in the IUE data is not due to extinction, but rather to the 
onset of broad emission from blended Fe n at 2300 A. 

For the neutral hydrogen intrinsic to NGC 3516 we find a best fit redshift of z = 0.0075±0.0005, 
Nhi = 3.5 x 10 17 cm -2 , and b = 10 km s _1 . This is blue-shifted by 400 ± 150 km s _1 relative to the 



2649 km s _1 systemic velocity of NGC 3516 measured using the stellar absorption lines (Vrtilek 
fe Carleton 1985] ). The opaque Lyman limit requires neutral hydrogen with a minimum column 
density Nhi > 2.2 x 10 17 cm -2 (90% confidence). Its sharpness limits b to less than 20 km s^ 1 at 
90% confidence. The effect of the assumed Doppler parameter on the shape of the intrinsic Lyman 
limit is shown with an enlarged view of the Lyman-limit region in NGC 3516 in Figure |. The 
general weakness of the Lyman absorption lines (only Lya and perhaps Ly/3 are detected) gives an 
upper limit of Nhi < 6.3 x 10 17 cm -2 at 90% confidence. 
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Fig. 3. — This expanded view of the short-wavelength end of the HUT spectrum shows the intrinsic 
Lyman limit and the O vi absorption more clearly. The position of the Galactic Lyman limit is 
shown to illustrate the intrinsic absorption. The solid, dotted, and dashed curves show the shape 
of the intrinsic Lyman absorption in NGC 3516 for assumed Doppler parameters of b = 10, 20, and 
50 km s _1 , respectively. As described in the text, at 90% confidence, b < 20 km s _1 . The heavy 
and light solid lines show our fit to the O vi emission and absorption. The light solid line shows 
the O vi emission profile superposed on the powerlaw continuum. The heavy solid line shows the 
contribution of the intrinsic O vi absorption to the fit. The light solid line in the vicinity of 1037 
A shows the contribution of intrinsic Ly/3 and Galactic C n A1037 to the fit. Galactic absorption 
features in the spectrum are denoted by a "(G)" following the line identification. 

Measured properties of the fitted emission lines are in Table [l]. The fitted absorption lines 
are summarized in Table |2[ All tabulated features, except for Ly/3, have a statistical significance 
exceeding 3a. The error matrix of the fit is used to derive the la statistical errors shown in the 
table. Systematic errors of up to ~5% are likely to be present in the fluxes, and our wavelength scale 
has a limiting accuracy of ~ 60 km s _1 . For ease in fitting, some poorly constrained parameters 
were linked to others with better determined values. For example, the velocity offset and FWHM 
of the C in A977 and the N in A991 lines were linked to share common values, and the FWHM of 
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He II A1085 was set equal to that of He n A1640. Any entries in the tables with identical velocity 
offsets or FWHM's were linked in a similar manner. Since the intrinsic O vi absorption is difficult 
to see in Figure [I], an enlarged view illustrating our fits to this region is shown in Figure ||. The 
O vi absorption doublet straddles the peak of the O vi emission line, and the sharpness of the line 
peak is largely defined by the flux escaping at wavelengths between the doublets. Galactic C II, 
Galactic Ly/3, and intrinsic Ly/3 absorption all contribute to the deficit in flux on the blue wing of 
the O vi emission. 



Table 1. Emission Lines in the HUT Spectrum of NGC 3516 



Line 


Avac 


Flux a 


Velocity 13 


FWHM 




(A) 




(km s^ 1 ) 


(km s _1 ) 


P TTT 


Q77 09 


U.l I u.o 




O-Liyu HZ 'iuy 


N III 


990.83 


2.3 ±0.3 


-235 ± 249 


3190 ± 469 


VI 


1033.83 


15.9 ± 1.1 


13 ± 63 


3516 ± 375 


He II 


1085.15 


2.9 ±0.6 


-2012 ± 880 


10041 ± 1694 


Lya 


1215.67 


5.6 ±0.9 


721 ± 33 


808 ± 98 


Lya 


1215.67 


84.3 ±3.1 


-345 ± 160 


4556 ± 230 


Lya 


1215.67 


89.9 ±3.2 






N V 


1240.15 


5.2 ± 1.0 


1443 ± 339 


5660 ± 181 


I 


1304.35 


1.3 ±0.3 


-1152 ± 573 


5660 ± 181 


C II 


1335.30 


1.9 ±0.3 


285 ± 103 


5660 ± 181 


Si IV 


1393.76 


6.1 ±0.2 


285 ± 103 


5660 ± 181 


Si IV 


1402.77 


3.1 ±0.1 


285 ± 103 


5660 ± 181 


N IV] 


1486.50 


1.1 ±0.2 


-49 ± 242 


2152 ± 498 


C IV 


1549.05 


58.4 ±0.6 


-191 ± 65 


3705 ± 117 


He II 


1640.50 


2.4 ±0.5 


-280 ± 113 


1869 ± 285 


He II 


1640.50 


5.7 ±0.9 


-280 ± 113 


10041 ± 1694 


He II 


1640.50 


8.2 ± 1.0 






N III] 


1750.51 


2.2 ±0.4 


-1084 ± 238 


2870 ± 476 


a Flux in 10" 


-13 —2 - 

erg cm s 


1 corrected for E(B 


- V) = 0.06. 




b Velocity relative to a systemic redshift of cz = 


2649 km s" 1 (Vrtilek & 


1 Carleton 1985). 
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3. The Broad Emission Lines 

The broad emission lines in NGC 3516 as observed with HUT show subtle differences compared 
to other Seyfert l's and low redshift AGN. The optical to X-ray spectral index for NGC 3516, 
a ox = 1.20, is typical of other Seyfert l's flKriss fc Canizares 1985] ), but the flux ratio of O vi 
A1034 to Lya is only 0.18. This lies below the correlation of O Vi/Lya with a ox discussed by 
Zheng, Kriss, & Davidsen (1995), but is comparable to the mean value of 0.17 seen in high redshift 
quasars ( [Laor et al. 1994| ). As in the case of Fairall 9, however, the lack of a soft X-ray excess in 



Table 2. Absorption Lines in the HUT Spectrum of NGC 3516 



Line 


/v vac 


EW 


Velocity a 


FWHM 




(A) 


(A) 


(km s 


x ) 


(km s _ ) 


N I 


953.77 


1.3 ±0.4 


-2944 ± 


146 


804 ± 274 


N I 


964.24 


1.6 ±0.5 


-2392 ± 


115 


804 ± 274 


C III 


977.02 


< 1.0 







800 b 


Si II 


989.87 


0.8 ±0.3 


-2704 ± 


107 


804 ± 274 


Ly/3 


1025.72 


2.1 ±0.4 


-2649 ± 


156 


880 b 


Ly/3 


1025.72 


0.2 ±0.4 


-400 ± 


150 


880 b 


C II 


1036.90 


1.2 ±0.3 


-2798 ± 


109 


1076 ± 146 


VI 


1031.93 


1.6 ±0.4 


-161 ± 


81 


1076 ± 146 


VI 


1037.62 


0.7 ±0.4 


-161 ± 


81 


1076 ± 146 


Si II 


1193.14 


2.2 ±0.4 


-2733 ± 


127 


1645 ± 255 


N I 


1200.16 


2.6 ±0.5 


-2269 ± 


122 


1645 ± 255 


Lya 


1215.67 


0.5 ±0.2 


-400 ± 


150 


770 b 


N V 


1238.82 


0.9 ±0.2 


63 ± 


53 


632 ± 116 


N V 


1242.80 


0.4 ±0.2 


63 ± 


53 


632 ± 116 


Si II 


1260.42 


1.3 ±0.3 


-2456 ± 


63 


831 ± 170 


C II 


1335.30 


1.6 ±0.3 


-2454 ± 


58 


789 ± 146 


Si IV 


1393.76 


0.8 ±0.4 


-204 ± 


75 


778 ± 585 


Si IV 


1402.77 


0.4 ±0.3 


-204 ± 


75 


778 ± 585 


Si II 


1527.17 


0.6 ±0.2 


-2495 ± 


47 


324 ± 106 


C IV 


1548.19 


1.9 ±0.1 


-15 ± 


22 


516 ± 37 


C IV 


1550.77 


1.0 ±0.1 


-15 ± 


22 


516 ± 37 



a Velocity relative to a systemic redshift of cz = 2649 km s 1 ( Vrtilek fc Carleton 1985| ). 
b FWHM fixed at the instrument resolution at this wavelength. 
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NGC 3516 may be a significant factor in producing lower O vi/Lya compared to other Seyfert l's 
flZheng et al. 1995|). 



The other noticeable difference is the relatively high strength of broad C in A977 and N in 
A991 emission in our spectrum. Although these lines have been seen in other low-redshift quasars 



(Laor et al. 1994; Laor et al. 1995), they are not detected in other HUT spectra of Seyfert l's. The 
ratios of these lines to each other and to O vi A1034 is also typical of that seen in the low redshift 
quasars observed with HST ([Laor et al. 1995| ). 

Reverberation mapping experiments have shown that the broad emission-line region (BELR) 



of AGN is highly stratified in both spatial distribution and in ionization parameter (e.g., Clavel et 



al. 1991| ; jKrolik et al. 1991| ; [Reichert et al. 1994] ; |Korista et al. 1995| ). Nevertheless, it is useful 



to match a single-zone photoionization model to the observed broad line emission in NGC 3516 to 
obtain a fiducial for comparison to other Seyferts. In addition, it provides a comparison to models 
of the warm absorbing medium seen in the X-ray spectrum and in the UV absorption lines that 
we discuss in §4. We use the photoionization code XSTAR ( Kallman &: Krolik 1993] ) to compute 



a grid of photoionization models varying the total column density N and the ionization parameter 
U = ni on /nj{, where nj on is the number density of ionizing photons between 13.6 eV and 13.6 
keV illuminating the cloud and nu is the density of hydrogen atoms. As we are constraining 
only high ionization lines that are insensitive to the density, we assume constant density clouds 
with Tin = 10 9 cm -3 and solar abundances. For the incident photoionizing continuum we use the 
extinction-corrected UV and absorption-corrected X-ray spectrum of NGC 3516 as described by 
Kriss et al. (1996). The UV power law was extrapolated to higher energies following f u ~ j/ -1 - 89 
with a break at 51 eV to the slope of the X-ray power law, f u ~ z;~ - 78 . To match the observed 
line ratios we varied the parameters until we achieved the closest fit to the strongest lines - 
O vi+Ly/3, Lya, Si iv+O iv], C iv, and He n A1640. The closest match is for an ionization 
parameter U = 0.035. The resulting line ratios from this model are compared to the observed 
values and their error bars in Table ||. Choosing U = 0.035 is mainly a balance between the O vi 
emission and the C iv line. Producing sufficient O vi requires U = 0.045, but then the C iv/Lya 
ratio becomes too high. U = 0.030 fits the C iv/Lya intensity ratio well, but the O vi emission is 
then a factor of 2 less than observed. As the lines we are matching originate in the high ionization 
illuminated faces of the BELR clouds, their intensities are rather insensitive to changes in the 
column density above log N = 21.6 in our models. 

The best-fit ionization parameter U = 0.035 is a typical value for single-zone models of AGN 
broad-line regions, but one can see from Table [3| that this model is insufficient to explain all the 
observed line ratios. The most noticeable differences are in the relatively high strengths observed 
for C in A977, N in A991, and He n A1640. Given the compromise we made between matching 
O vi and C iv, one can envision that a higher ionization zone producing relatively more O VI and 
He II A1640 could account for part of these differences. An additional population of higher density 
clouds may be required to produce the enhanced C in A977, which becomes a more important 
coolant as other lines become optically thick ( Netzer 1990j ). The strong N in A991 is more of a 
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puzzle. Most photoionization models predict it to be weaker than C in A977, yet in our spectrum 
and in the HST quasar spectra ( [Laor et al. 1995| ) it is stronger. It is possible that fluorescent 



mechanisms could enhance the strength of this transition under favorable circumstances (Ferguson 
Ferland, fc Pradhan 1995| ). 

Although the broad lines in NGC 3516 seem to be mostly produced in a region with an 
ionization parameter typical of the BELR in other AGN, U = 0.035 is an order of magnitude lower 
than the ionization parameter required for the warm absorbing gas detected in the ASCA X-ray 
spectrum ( Kriss et al. 1996 ). This does not rule out an origin in the BELR for the X-ray warm 
absorbers, but these absorbers must be physically distinct from the clouds producing the bulk of 
the broad-line emission. 



4. The Complex Absorbing Medium in NGC 3516 

Intrinsic UV absorption lines are prominent in our spectrum of NGC 3516, but they are near 
the weakest levels seen in NGC 3516, comparable to their appearance in the 1993 IUE monitoring 



Table 3. Observed and Modeled Emission Line Ratios in NGC 3516 



Feature 


^/^A1216 a 


Model b 


C in A977 


0.008 ± 0.003 


0.002 


N in A991 


0.026 ± 0.004 


0.000 


O vi+Ly/J 


0.177 ±0.014 


0.123 


He ii A1085 


0.032 ± 0.007 


0.009 


N v A1240 


0.058 ±0.011 


0.064 


O i A1304 


0.014 ± 0.003 


0.000 


C ii A1335 


0.021 ± 0.003 


0.007 


Si iv+O iv A1400 


0.102 ±0.005 


0.102 


N iv] A 1486 


0.012 ±0.002 


0.074 


C iv A1549 


0.650 ± 0.024 


0.778 


He ii A1640 


0.091 ±0.012 


0.037 


N in] A1750 


0.024 ± 0.004 


0.012 



a Line ratios are relative to the major contributors to the A1216 blend — H Lya, He n 
A1216, and O v A1218. 

Line ratio computed for a model with U = 0.035, n H = 10 9 cm , and ./V = 10 cm . 
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campaign flKoratkar et al. 1996 ), The frequently observed blue-shifted absorption trough 



5 hull, k. Begelman 1987| ; [Walter et al. 1990| ; |Kolman et al. 1993| ) is not present, and the C iv 



feature in our spectrum probably corresponds to the narrow component in the models of Walter et 
al. (1990). Although the C iv absorption line in the HUT spectrum is narrow, it is resolved. After 
correcting for instrumental broadening by subtracting a 2.0 A Gaussian in quadrature (Kruk et al 



1995 ), we find an intrinsic width for the line of 350 km s _1 . This corresponds to b = 210 km s _1 , 
which is consistent with the observed equivalent widths (EWs) and optically thin doublet ratio of 
the deblended C IV lines. The other high ionization resonance lines also have optically thin doublet 
ratios consistent with b ~ 200 km s _1 . In contrast, the neutral hydrogen absorption at the Lyman 
limit requires b < 20 km s _1 . This, and the significantly different velocity of the neutral hydrogen 
absorption, are indications that the absorbing medium contains multiple zones. 

The simultaneous X-ray observations of NGC 3516 discussed in the preceding paper flKriss et 



al. 1996) also hint that the absorbing medium may be complex. The ASCA observations show 



O vii and O viii edges characteristic of warm absorbers that have been seen in other AGN (Turner 



et al. 1993| ; pVlathur et al. 1994] ; |Fabian et al. 19944 |Fabian et al. 1994b| ). Photoionized warm 



absorber models similar to those discussed by Krolik & Kriss (1995) can fit the data, but the fit is 
significantly improved if two absorbing zones with differing ionization parameter and total column 
density are used. These warm absorber models use the same incident photoionizing continuum 
described in the previous section. We also considered alternative ionizing continua to test the 
sensitivity of our results to the assumed incident spectrum. We used both an extremely hard 
spectrum with f u ~ j,— 0.78 f rom 2500 A through the UV and X-ray up to 100 keV, and also the 
spectral shape of NGC 5548 as used by Mathur et al. (1995), which contains a soft X-ray excess. 
The models are computed for constant density (n# = 10 3 cm -3 ) clouds in thermal equilibrium. 
(For densities < 10 11 cm" 3 there are no density-dependent effects in our calculations or results.) 
The models are described by their ionization parameter U and their total column density N. 



Single-zone photoionization models of the absorbing gas in 3C 351 (Mathur et al. 1994) and 



in NGC 5548 (Mathur, Wilkes, & Elvis 1995) can produce both the X-ray and UV absorption seen 
in these two AGN. We now test whether this is true for NGC 3516. The single warm absorber fit 
to the ASCA spectrum of NGC 3516 has U = 0.48 and logiV = 22.02. Given the column densities 
predicted by this model for the UV ions, we can place the observed EWs on curves of growth 
(Figure Q). If the model self-consistently accounts for both the X-ray and UV absorbers, all the 
plotted EWs should lie on a single curve of growth. As shown in the top panel of Figure |4], they do 
not. The predicted column densities of neutral hydrogen and Si iv are orders of magnitude below 
what we observe in the HUT spectrum, and a single Doppler parameter does not give a good match 
to the remaining UV lines. Choosing b = 70 km s _1 fits the N v and O vi doublets within the 
errors, but this underpredicts the strength of the C iv doublet by factors of several. This situation 
is similar to the difficulties noted by Kriss et al. (1995) in trying to account for the low-ionization 
UV absorption lines in NGC 4151 with a single warm absorber model. 



The warm absorber models computed with the alternative ionizing continua give similar results, 



-12- 




14 16 18 20 

log N,fA (cm 2 A) 



Fig. 4. — Top Panel: The observed EWs of the UV absorption lines in NGC 3516 are plotted on 
curves of growth using column densities predicted by the single warm absorber fit to the ASCA 
X-ray spectrum (Kriss et al. 1996). This model has U = 0.48 and a total column density of 
10 22 - 02 cm~ 2 . Points are plotted at a horizontal position determined by the column density for 
the given ion in the model with a vertical coordinate determined by the observed EW for the 
corresponding absorption line. The vertical error bars are from Table 2, and the horizontal error 
bars are the range in column density allowed by the uncertainty in the fit to the ASCA spectrum. 
The predicted column of Si iv in this model is too low to appear on the plot. The thin solid lines 
show predicted EWs as a function of column density for Voigt profiles with Doppler parameters 
of 6 = 10, 20, 50, 100, 200, and 500 km s" 1 . A model that fits the data would have all points 
lying on one of these curves. This model cannot simultaneously match both the UV and the X-ray 
absorption. Center Panel: The observed EWs are plotted for column densities predicted by a model 
with U = 0.025 and a total column density of 10 19 ' 2 cm -2 . The heavy solid curve at b = 200 km s _1 
gives a good match to the observed EWs of the C iv, N v, and O vi doublets, but the predicted 
EW of Si iv is far below the observed value. Bottom Panel: The observed EWs are plotted for 
column densities predicted by a model with U = 0.0013 and N = 10 19 cm~ 2 . b = 20 km s _1 can 
match the observed EWs of Si iv and Ly/3 and satisfy the upper limit on the EW of C in A977. 

with the predicted columns of the UV ions in each model within tens of percent of each other. This 
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insensitivity to the precise shape of the ionizing spectrum is not surprising. For broad ionizing 
continua such as power laws or broken power laws, even the earliest photoionization calculations 
(e.g., [Tarter, Tucker, fc Salpeter 1969| ) showed that models were much more sensitive to ionization 
parameter than to spectral shape. In Table |I| we compare the columns inferred from the UV 
absorption lines assuming they are optically thin, as indicated by their doublet ratios, to the 
column densities predicted by the various models. The similarity of the models to each other is 
apparent, as is the large disagreement with the observations, particularly for Si iv. 

The two-component warm absorber fit favored by the ASCA data also fails as an adequate 
description of the UV absorption. The lower ionization zone with U = 0.32 and logiV = 21.84 
still falls short in matching the Si iv and H I by orders of magnitude, and once again a single 
Doppler parameter fails to match the observations of the remaining UV lines. The only UV- 
absorbing species with a significant column at any Doppler parameter in the higher ionization zone 
(U = 1.66, logiV = 22.15) is O VI. 

Producing the observed EWs of the UV lines requires regions of lower ionization and lower 
column density than those producing the X-ray absorption. For U = 0.025, logiV = 19.2 cm -2 , 
and b = 200 km s -1 we find a good match to the observed EWs of the C iv, N v, and O vi 



Table 4. Observed and Predicted Absorption Columns in NGC 3516 



Ion 


N obs a 


NModell h 


^ModeU 


NModel3 d 




(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


H I 


3.5 x 10 17 


5.0 x 10 16 


1.7 X 10 17 


1.8 x 10 16 


C iv 


4.7 x 10 14 


4.7 x 10 14 


6.0 x 10 14 


5.4 x 10 14 


N v 


4.1 x 10 14 


2.1 x 10 15 


2.5 x 10 15 


3.2 x 10 15 


O vi 


1.3 x 10 15 


9.8 x 10 16 


1.1 x 10 17 


1.4 x 10 17 


Si iv 


1.1 x 10 14 


< 1 x 10 10 


< 1 x 10 10 


< 1 x 10 10 



a Column density assuming the observed absorption lines are optically thin. The H I column is 
from the Lyman-limit fit. 

b Predicted column density for the single-zone warm absorber model with the NGC 3516 ionizing 
spectrum, U = 0.48 and log N = 22.02. 

c Predicted column density for the warm absorber model with the j/~ - 78 ionizing spectrum, U = 
0.16 and log N = 22.02. 

d Predicted column density for the warm absorber model with the NGC 5548 ionizing spectrum, 
U = 1.47 and log N = 22.02. 
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doublets as shown in the center panel of Figure ||. The predicted EW of the Ly/3 line is ~ 2a higher 
than observed, but enhancing the overall metal abundance by a factor of a few can correct this. 
However, the predicted EW of Si iv in this model is orders of magnitude below the observed value. 
Similar problems are encountered in matching the inferred column densities of Si iv absorption in 
models of broad-absorption line (BAL) QSOs ( Weymann, Turnshek, fc Christiansen 1985| ). Ad hoc 



adjustments to the ionizing continuum can alleviate the disagreement to some extent (|Weymann 



|Turnshek, &: Christiansen 1985| ), but even the most extreme changes do not bring the predictions 
within the bounds of the observations, again showing the relative insensitivity of photoionization 
models to the shape of the ionizing spectrum. Producing enough Si iv absorption requires a zone 
of even lower ionization. Choosing U = 0.0013 and logiV = 19.0 matches the Si iv EWs without 
contributing much to the columns of the higher ionization lines. This model is illustrated in the 
lower panel of Figure ||. Limits on the EWs of lower ionization species such as Ly/3 and C in 
A977, however, require b < 50 km s _1 . This lower ionization zone with low Doppler parameter 
is potentially the origin of the optically thick Lyman limit system, but the predicted column of 
neutral hydrogen in this model is only Njji = 7.7 x 10 16 cm -2 , a factor of 4 below the observed 
value. 

The multiple zones considered above span a range of 10 3 in both column density and ionization 
parameter, and they form a nearly orthogonal set that makes it possible to consistently model the 
UV and X-ray absorption with separate regions of gas in very different physical states. For Doppler 
parameters at the sound speed of the two high ionization zones producing the X-ray absorption 
(30 km s~ x in the lower ionization zone, 70 km s _1 in the higher), the predicted contribution of 
the X-ray absorbing gas to the EWs of the C iv, N v, and O vi lines is only 20-30% of their 
observed values. The UV zones producing the C iv, N v, O vi and Si iv lines have total columns 
with negligible impact on the X-ray opacity. Since the X-ray absorbing gas is predicted to have a 
significantly lower Doppler parameter than the gas producing the bulk of the UV absorption lines, 
high resolution spectra in the UV should be able to identify kinematic components associated with 
the different zones. 

Several previous authors have suggested that the UV absorption lines and the X-ray absorption 



seen in some AGN may have a common origin in BELR clouds ( Ferland fc Mushotzky 1982 ; Re 



ichert, Mushotzky, fc Holt 1986| ). However, Mathur et al. (1994) show that the physical conditions 



of the UV and X-ray absorber in 3C 351 are different from those of the BELR clouds, with the 
X-ray warm absorber requiring a much higher ionization parameter than the clouds producing the 
broad emission lines. For NGC 3516, Voit, Shull, & Begelman (1987) argued that the observed 
absorption profile (the C iv absorption in the 1980's was much broader and stronger than observed 
here) required 10-30 BELR clouds along the line of sight and that photoionized C iv could not be 
present in the outer clouds due to shielding by the inner ones. The lack of observed Mg n absorption 
was another difficulty for models involving BELR clouds. They concluded that the absorption must 
arise in optically thin, outflowing material with a density exceeding 10 5 cm -3 , given the observed 
variability, and that it could not be associated with typical BELR clouds. 
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Our conclusions are similar, even though the absorption is now much weaker. The ionization 
parameter we derive for the broad emission-line clouds (U ~ 0.035) is an order of magnitude lower 
than that required by the X-ray absorbing gas, making it highly unlikely that the X-ray absorption 
occurs in BELR clouds. U ~ 0.035, however, is similar to the ionization parameter we find for the 
zone producing the C iv, N v, and O vi absorption lines (U ~ 0.025). Although these absorption 
lines and the broad emission lines share comparable ionization parameters, the total columns of 
the BELR clouds are ~ 10 3 x higher, as the broad emission line clouds are optically thick in 
these transitions, but the absorbing clouds are optically thin. This column density constraint is 
a poor one, however, because the BELR is highly stratified, and there may well be a population 
of clouds that both make the absorption lines and contribute to the broad-line emission. The 
velocity distribution of such clouds would have to be rather peculiar, however, to project a velocity 
width of only 200 km s _1 along the line of sight from an ensemble with an emission line width 
of ~ 4000 km s _1 . We conclude that we still cannot unambiguously establish the location of the 
absorbing gas, but it must be physically distinct from the broad emission line clouds. 

Although we have described the UV and X-ray absorbing zones as discrete entities, it is quite 
likely that there is a broad, possibly continuous, distribution of parameters such as one might find 
in an outflowing wind, either from the accretion disk or from the surface of the obscuring torus. 
As discussed in the introduction, NGC 3516 is one of the rare Seyfert 1 galaxies that shows an 
extended, biconical NLR (see Golev et al. 1995 for the best images). Among Seyfert l's, only 
NGC 4151 has such an extensive NLR with a biconical morphology (|Pogge 1989| ; Evans et al. 



1994 ; Schmitt fc Kinney 1996 ). The biconical morphology in NGC 3516 implies that our line of 
sight passes close to the surface of the obscuring torus as suggested for NGC 4151 by Evans et al. 
(1993). The opaque Lyman limit is another characteristic that NGC 3516 shares with NGC 4151. 
These features are unusual in combination for a Seyfert 1, and their presence in yet another galaxy 
strengthens the case for collimation of the ionizing radiation by the absorbing gas. 

The complexity of the UV and X-ray absorption that both NGC 3516 and NGC 4151 exhibit 
suggest that the inclination of the source relative to the observer may lead to the differences between 
these two objects and the apparently simpler cases considered by Mathur et al. (1994,1995). The 
simpler warm absorbers with only high-ionization UV absorption lines may be viewed at higher 
inclination, further away from the denser medium near the torus. Such geometrical differences may 
ultimately help us to understand the location and origin of the warm absorbing gas in AGN. 
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